Divergent spike mutations impact the activation of the fusion core in Delta and Omicronvariants of SARS-CoV-2 (preprint)

SARS-CoV-2 infects host cells by binding the receptor-binding domain (RBD) of its spike protein to the receptor, ACE2. A subset of highly effective spike mutations plays critical roles in altering the conformational dynamics of spike protein. Here, we use molecular dynamics simulations to investigate how spike mutations affect the conformational dynamics of spike/ACE2 complex in the D614G, Delta (B.1.617.2) and Omicron (B.1.1.529) SARS-CoV-2 variants. We observe that the increased positive-charged mutations in the Omicron spike amplify its structural rigidity and reduce its structural flexibility. The mutations (P681R in Delta and P681H in Omicron) at the S1/S2 junction facilitate S1/S2 cleavage and aid the activation of the fusion core. We report that high structural flexibility in Delta lowers the barrier for the activation of the S2 core; however, high structural rigidity in Omicron enhances the barrier for the same. Our results also explain why Omicron requires the presence of a higher number of ACE2 to activate its fusion core than Delta.

The SARS-CoV-2 nucleoprotein associates with anionic lipid membranes (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a lipid-enveloped virus that acquires its lipid bilayer from the host cell it infects. SARS-CoV-2 can spread from cell to cell or from patient to patient by undergoing assembly and budding to form new virions. The assembly and budding of SARS-CoV-2 is mediated by several structural proteins known as envelope (E), membrane (M), nucleoprotein (N) and spike (S), which can form virus-like particles (VLPs) when co-expressed in mammalian cells. Assembly and budding of SARS-CoV-2 from the host ER-Golgi intermediate compartment is a critical step in the virus acquiring its lipid bilayer. To date, little information is available on how SARS-CoV-2 assembles and forms new viral particles from host membranes. In this study, we find the N protein can strongly associate with anionic lipids including phosphoinositides and phosphatidylserine. Moreover, lipid binding is shown to occur in the N protein C-terminal domain, which is supported by extensive in silico analysis. Anionic lipid binding occurs for both the free and N oligomeric forms suggesting N can associate with membranes in the nucleocapsid form. Herein we present a lipid-dependent model based on in vitro, cellular and in silico data for the recruitment of N to M assembly sites in the lifecycle of SARS-CoV-2.

Structure of SARS-CoV-2 M protein in lipid nanodiscs (preprint)

SARS-CoV-2 encodes four structural proteins incorporated into virions, spike (S), envelope (E), nucleocapsid (N), and membrane (M). M plays an essential role in viral assembly by organizing other structural proteins through physical interactions and directing them to sites of viral budding. As the most abundant protein in the viral envelope and a target of patient antibodies, M is a compelling target for vaccines and therapeutics. Still, the structure of M and molecular basis for its role in virion formation are unknown. Here, we present the cryo-EM structure of SARS-CoV-2 M in lipid nanodiscs to 3.5 [A] resolution. M forms a 50 kDa homodimer that is structurally related to the SARS-CoV-2 ORF3a viroporin, suggesting a shared ancestral origin. Structural comparisons reveal how intersubunit gaps create a small, enclosed pocket in M and large open cavity in ORF3a, consistent with a structural role and ion channel activity, respectively. M displays a strikingly electropositive cytosolic surface that may be important for interactions with N, S, and viral RNA. Molecular dynamics simulations show a high degree of structural rigidity and support a role for M homodimers in scaffolding viral assembly. Together, these results provide insight into roles for M in coronavirus assembly and structure.

Cooperative multivalent receptor binding promotes exposure of the SARS-CoV-2 fusion machinery core (preprint)

The molecular events that permit the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to bind, fuse, and enter cells are important to understand for both fundamental and therapeutic reasons. Spike proteins consist of S1 and S2 domains, which recognize angiotensin-converting enzyme 2 (ACE2) receptors and contain the viral fusion machinery, respectively. Ostensibly, the binding of spike trimers to ACE2 receptors promotes the preparation of the fusion machinery by dissociation of the S1 domains. We report the development and use of coarse-grained models and simulations to investigate the dynamical mechanisms involved in viral binding and exposure of the S2 trimeric core. We show that spike trimers cooperatively bind to multiple ACE2 dimers. The multivalent interaction cyclically and processively induces S1 dissociation, thereby exposing the S2 core containing the fusion machinery. Our simulations thus reveal an important concerted interaction between spike trimers and ACE2 dimers that primes the virus for membrane fusion and entry.

Molecular interactions of the M and E integral membrane proteins of SARS-CoV-2 (preprint)

Specific lipid-protein interactions are key for cellular processes, and even more so for the replication of pathogens. The COVID-19 pandemic has drastically changed our lives and cause the death of nearly three million people worldwide, as of this writing. SARS-CoV-2 is the virus that causes the disease and has been at the center of scientific research over the past year. Most of the research on the virus is focused on key players during its initial attack and entry into the cellular host; namely the S protein, its glycan shield, and its interactions with the ACE2 receptors of human cells. As cases continue to raise around the globe, and new mutants are identified, there is an urgent need to understand the mechanisms of this virus during different stages of its life cycle. Here, we consider two integral membrane proteins of SARS-CoV-2 known to be important for viral assembly and infectivity. We have used microsecond-long all-atom molecular dynamics to examine the lipid-protein and protein-protein interactions of the membrane (M) and envelope (E) structural proteins of SARS-CoV-2 in a complex membrane model. We contrast the two proposed protein complexes for each of these proteins, and quantify their effect on their local lipid environment. This ongoing work also aims to provide molecular-level understanding of the mechanisms of action of this virus to possibly aid in the design of novel treatments.

A Multiscale Coarse-grained Model of the SARS-CoV-2 Virion (preprint)

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations, however, are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and on-going development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Structural data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data becomes publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion. Significance StatementThis study reports the construction of a molecular model for the SARS-CoV-2 virion and details our multiscale approach towards model refinement. The resulting model and methods can be applied to and enable the simulation of SARS-CoV-2 virions.